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Casimir energy changes are investigated for geometries obtained by small but arbitrary deforma- 
tions of a given geometry for which the vacuum energy is aheady known for the massless scalar 
field. As a specific case, deformation of a spherical shell is studied. From the deformation of the 
sphere we show that the Casimir energy is a decreasing function of the surface to volume ratio. The 

- - - decreasing rate is higher for less smooth deformations. 
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I. INTRODUCTION 



Casimir energies are known for several cavities in several spatial dimensions for electromagnetic or massless scalar 

fields. Exact Casimir energy calculations are available for rectangular prisms 1], for spherical shell 0, for cylindrical 

region i3|, for a pyramidal cavity and for a conical cavity [J]. All these geometries have definite boundary wall shapes. 

—^ , For example the prisms are all with right angular wedges, the sphere is the perfect one; and, the pyramid and the 

flN^' cone are of very special types. The obvious reason of these restrictions is the fact that these are the regions for which 

(~| ] one can calculate the exact field modes with the required boundary conditions. Limited number of examples ( all 

'— ' . with rigid walls ) do of course not give much hint about the dependence of the Casimir energy on the shapes of the 

^^ ; regions. 

> 

/-o There are to our knowledge two approaches in dealing with rather arbitrary geometries. First one is the 

QQ proximity force approximation [5j]. It is applicable to two body systems which are close to each other. It employs 
C^ the parallel plate modes in every cylindrical region of infinitesimal base between the bodies and then integrates 
\l \ over these regions [6|. Second approach is the one called multiple scattering expansion. It formulates the vacuum 
^+ ■ energy for the electromagnetic field in terms of the successive scattering from the conducting boundaries |7| . In this 
C^ ' approach between two successive scatterings free Green function is employed. The method enables one to investigate 
00 . the connections between the divergencies and the geometrical details of the boundaries. Employment of the free 
^^ ' Green functions however from one scattering to the next one is not of much practical value in compact space regions: 
^ , Specially for more curved boundaries one may need to consider large number of scatterings to approximate the exact 
Green function. On the other hand it may be more reasonable to employ the exact Green function ( if it is already 
known ) of the compact region ( instead of the free one ) if the region under consideration is sufficiently close to the 
jd original region. Variations around the exactly solvable geometries may give some hint about the dependence of the 
vacuum energies on the shape of the boundaries. 

In present work we try to investigate the shape dependence of the Casimir energy for massless scalar field by 
calculating the effect of the small but arbitrary deformations of a given geometry for which we already know the 
vacuum energies. We use Green function to formulate the perturbation theory around the exact solution of the 
region with boundary 5'. Suppose Gf is the exact Green function for the massless scalar field confined to the region 
in 5', and (3 is the small deformation of S. Converting a boundary problem into an integral equation we arrive at the 
perturbation series 

Gf = Gf + /3Gi^>/?2Gf + ••• 
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where Glf is the correction to the original Green function Gf resulting from the j times reflections from the deformed 
boundary S. One reflection gives information about the size of the new boundary. To get information about the 
shape dependence we need to take into account at least two reflections. 

Having in hand the Green function we can construct the zeta function which is useful tool in Casimir energy 
calculation. The zeta function of the system can be expressed in terms of the heat kernel coefficients which are 
functionals of the geometrical invariants of the boundaries Q. For arbitrary geometries these coefficients are too 
complicated and therefore are not available for practical purposes |9|. From this point of view we think that the 
perturbation around a known geometry may be an effective approach. 

For massless fields which are the only fields for which the Casimir energy is meaningful, there is no unique 
way of getting rid of the infinities if the heat kernel coefficient 02 of the zeta function expansion is not zero [l^l ■ The 
situation can be improved if one considers the whole space for when one sums the zeta functions of the in and out 
regions 02 coefficients cancel each other. Of course for such cancelation the boundaries should be free of sharp corners. 

In the coming section we briefly review the zeta function approach to the vacuum energy calculations and 
the regularization scheme which we employ in the our work. 

In Section III we present the general formulation of the Casimir energy contribution of small deformations of 
the boundaries. 

In Section IV deformation of the sphere is discussed. 

In Section V we analyze the dependence of the energy on the shape of the boundary. 

Details of the involved calculations are given in the Appendices. 

II. A BRIEF REVIEW OF THE ZETA FUNCTION METHOD 

Formally the calculation of the Casimir energy is reduced to a treatment of a sum over all one particle energy 
eigenvalues 



^ = ^EV^ (1) 
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This sum is divergent and regularization is needed. For the scalar field confined in a compact three dimensional region 
with the Dirichlet boundary condition the zeta function 



c(z) = ^i?r (2) 



n=l 



is well defined by Weyl theorem for Rez > 3/2 Q. By the analytic continuation it is possible to define this function 
on the whole complex plane. This may be done by using the representation 



C(z) = — - / dte-'K{t) (3) 

r(z) Jo 

with the heat-kernel 

00 
i^(t) = ^e-*^". (4) 

n=l 

For i ^ 00 the integral is well behaved. Possible poles arise from t ^ behavior of the heat kernel Q 

00 
K{t)^ J2 a^t"-'/'. (5) 

n=0,l/2,l,... 



Splitting the integral as Jg dt + J-^ dt we arrive at 

Res{C{z)T{z)) U=3/2-„= a„ (6) 

where a„ are heat kernel coefficients which depent on the geometry of boundary which confines the scalar field. 
When the coefficient 02 is nonzero the value of the zeta function at z = —1/2 which defines the Casiniir energy 
becomes infinite. For niassless scalar field there is no renormalization condition to get rid of this pole in a unique 
way [10| . One needs additional considerations to apply the zeta function method to the calculation of the vacuum 
energy. The extrinsic curvature of the sphere will have opposite sign when viewed from inside or outside. 02 heat 
-kernel coefficient depends on an odd power of extrinsic curvature. Two of them add to cancel each other when we 
approach the surface both from in and out region of the ball. This does not hold only for the spherical shell but 
is a general property for boundaries of an arbitrary shape. This cancelation of poles occurs only for infinitely thin 
boundaries. Once a finite thickness is introduced the absolute value of the extrinsic curvature at the inner and outer 
side of the boundary is different and divergencies do not cancel each other. 

In the present work we try to investigate the shape dependence of the Casimir energy for the massless scalar 
field by calculating the effect of the small but arbitrary smooth deformations of the boundary of given regions for 
which we already know the vacuum energies. We restrict our attention to the deformations of the spherical shell. 
Using the scattering theory in the spherical coordinates one arrives at the zeta functions inside 

"" 1=0 ■'^ ^ 

and outside the ball [3, ^^ 

CXS PCX) 1 

C<,„,(z) = ^^(2Z + l)y dojuj-'^-^HJ+'/'K^+y.iuj)) (8) 

The zeta function in the whole space 



/=o -^0 



az) = l{C^niz) + Cout{z)) (9) 

or 

^ t^ Jo dw ' 

is well defined at z = — 1/2. To find this value we use the uniform asymptotic expansions for Bessel functions |ll| 

^«-'-y^(T^<' + B->'^) (11) 



k=l 



and 



fc=i 
where 



u-^)- ^ /' (1+f:^) (12) 



t=^=L=, rj=Vl + x^+ln- ^== (13) 

Vl + x^ 1 + Vl + a? 



and coefficients Uk(t) satisfy the recurrence relation 



uu+i{t)^^t\\~e)u'^(t) + tj dr(l-5r2)ufe(r) (14) 



with the initial condition Mo(i) = 1- The uniform asymptotic expansions imply 



ln/.(^x)=f;^ (15) 



-1 



and 



lnif.M=f;(-)'=^ (16) 



where the first four terms of X„(i)s are 



X_i = i+ln- * 



t 1 + t 

t 5t3 

The zeta function in the whole space becomes 

C(-l/2) = - V Co (2m - 2) / da;X2™(t) (18) 

where 



oo 



■^ (1/2 + n) 

n=0 ^ ' ' ^ 

is the Riemann zeta function which vanishes for z = and z = — 2 [13|. Therefore the expansion (18) converges. The 
vacuum energy for the massless scalar field in the spherical shell of radius R is [iJl 

^sp/. = ^C(-l/2)^|, a ^0,003. (20) 

III. CONTRIBUTION OF SMALL BOUNDARY DEFORMATIONS TO THE VACUUM ENERGY 

In this section we use the Green function representation of the zeta function for the massless scalar field in the 
three dimensional space vanishing on a surface S 

^s^^^^M^ r^^^-2.+i f d3fGf(a-,f) (21) 

'^ Jo Jr3 

where 

rS(^,-t\-S G^""^(f,y), X,y€n,n (r,r,. 

^-^""^y' - ^ G°J'Hx,y), x,y€n,ut ^^^> 

is the Green function in R^ satisfying the boundary problem 

{-A + u;')G^Jx,y}^S{x-y}, Gf(x,y) = 0, x e S. (23) 

Here G™^ ( G°"*'^ ) is the Green function in the in-region flin ( the out-region flout ) • To have a well defined integral 
over R^ in (I2ip one considers a ball of radius L and then let it go to infinite. Divergent terms in powers of L will 
correspond to the infinite vacuum oscillations of the free Minkowski space. The variation 

<5C(^) = C^(z)-C^(z) (24) 

resulting from the deformation of the boundary S is then given by 

g^^^^^^MT^ r ^^^-2z+i j d^.^sGUx,x) (25) 

'^ Jo Jfl3 



where 

5G^{x, y) = Gt{x, y) ~ Cf (f , y). (26) 

Due to ([23]) it satisfies the wave equation 

{-A + oj^)5G^{x,y) = Q (27) 

and the boundary condition 

5G^{x,y) = -Gl{x,y), xeS. (28) 

The above boundary problem is equivalent to the integral equation [14| 

Gf (x, y) = Gf (x, y) - / dS^^^^G^M, v), (29) 



97n(w) 
where 



^ -(-li (30) 



dm{v) dv 

is the derivation along the unit vector rn{i)) normal to the wall S* at a point v. In the parametric representation 
V = v{t), t = {t^.t"^) the integration measure on S is 

ds = ^/\Y\(fT (31) 

where | 5 | is the determinant of the induced metric 

dv dv 

9ab={-g^,-g^) (32) 

with (•, •) being the scalar product in the three dimensional space. 

The solution of the integral equation (P5|) up to the second order ( which we can also interpret as the second 
reflection ) is 

,.aG5(f,.^)^s,,,^,^ f ,, f ,,dG^A^dG^{v^^s,J 



Gtix, y) = Gtix, y) - i d-s^^Gtiv, y) + j^ d~s j ^ '~^^^^^^>'^ V^ (^3) 



The property 



d_ 



d'xG'jz, a;)Gf (x, P) = ~^-^Gf (z, z') (34) 



allows us to integrate explicitly the perturbation solution over the three dimensional spatial space to get 

2 ouj"^ J§ om(v) aw^ J§ Jg om(v) dm[vi) 

We consider deformations of the boundary S along the unit vector n{z) normal to the surface 5* at a point z: 

v = z-(3n{z)f{z) (36) 

where /3 is the dimensionless deformation parameter of the surface S. This deformation formula implies 

9ab = gab-l3{g^,-g^) + ia^b)) (37) 

where gab is the metric tensor on 5*. Using (5-\/| g \ — l/2-y/| g {g'^'^Sgab we arrive at the variation of the integration 
measure 

d~s = ds-Pg^\§,,^)ds. (38) 



which together with the Taylor expansion 



.5 ,^ ,-.^ _ a .,^ 9Gt {x, z) , /3%2 , ^ d^Gt (f , z) 



imphes 



Gii^.^ = -pf^^^^ + T^ (^^)^#pr + ■ ■ ■ ^''^ 



Aduj^Js dn^{z) duj^ Js Js dn{z)dn{z') 



Up to the second order in (3 the zeta function variation ([24]) becomes 

sin(7rz) r ^^..-2z+i, 



5C{z) = — ^^ / dLJUj-^^+'5G^. (41) 



IV. DEFORMATION OF THE SPHERICAL SHELL 



The in and out-Green function for the massless scalar field vanishing on the sphere of the radius R are 

1=0 

and 



G'j;{r,n;r'n') = ^ V(2/ + l)P;((n, n')) ' — ' '+i/-^V J — (+V2V ) i+i/-\ /_ ^^3) 

AirVrr' frf, Ii+i/2[^R) 



Gr(r,n;r'n') = ^ V (2/ + l)PK(n, n')) ^K^-^M-R)h^„2{-r) ^ /.+r/2(a^i^)/^.+i/2(c.0)i^z+i/.(^/) ^^3^ 

47rVrr' -^ Ki+i/2[ujR) 



where < r < r' < R, Pi{x) is the Legendre polynomial and n is the unit vector normal to the sphere ( see Appendix 

d ^ _d_ 

dn{z) dr ' 



B ). The derivative normal to the sphere is 0-7^ = ■§-■ We have the derivatives of the following type 



1=0 



where 



Ti+^{^) = i^ln(/,+ i(w)J^,+ i(c.)) (45) 

is the spectral function in the whole space. Inserting the above type terms in ([3D]), (|4ip becomes 

1=0 -^o 

+ ^^^^ E (' + ^)(^' + ^)A. rci..^i-^^(r,+ .(^)-r,+.(^))2 (46) 

Ar = j^^ J dn I dn'Pi{in,n'))Pi,{{n,n'))f{n)f{n'), (47) 

dfJ = d(j)ddsm6 is the integration measure on the sphere. The first term in ([46]) at z = —1/2 is proportional to the 
vacuum energy of the massless scalar field confined in the spherical region of the radius R. Therefore the variation of 
the vacuum energy is then 



where 



/O r,2 9 °° 1 1 f°" 

dn{f-f{n) + f--f{n)) + — ^(/+ )(r + -)A, / du;cj'{T,^.{u;)^T,^. 
"" "" "" i,i'=o •^° 



{u;)r (48) 



The expansion (jBSp and the addition formula (|B4|) imply 

,7 

(2J+1)2 



■■' I cJ 12 



A.==/3^E E ^-T^Ml (49) 



J=OA/=-J 



where iiT;!, are the Clebsch Gordon coefRcients (|B6P and //^ are the expansion coefhcients of the deformation function 
/ in the spherical harmonics: 



oo J 



/(^) = E E /m^/C")- (50) 

J=0 A/=-J 

Using (|49|) we may represent (|48l) as 

^^ = Esphi^ J dQfin) + ^J dnifin) + f{n)Hf{n))] (51) 

where iJ is an energy operator 

HYi,{n)^H{J)Yij{n) (52) 

with e-values 

1 _^^'^_ 1 "^ /'OO 

H{J) = —Y.{l + -)Y.^NG■l,{^l) dLuu\Ti^.{oj)~Ti^^^.{Lo)f. (53) 

Here we used the formulae (IB6p and (|B12p for the Clebsch Gordon coefficients. The evaluation of H{J) which is 
quite involved, is given (together with all auxilary formulas ) in the Appendices. The energy operator H{J) has the 
following expansion 

H{J)^a^J^+a2J^+a^J+^^ + -^^^\-^ + --- . (54) 

J + 2 {J + 2) 

Terms with coefficients a_„, n = —1, —2, —3, ... in the above expansion give negligible contributions to the Casimir 
energy compared to the first three ones 

iJ(J) =-3,03 J3- 3,37 J2- 0,52 J. (55) 

The Casimir energy in the cavity obtained by the small but arbitrary deformations of the spherical region of radius 
R is, therefore given by 



R^ AttR / ■' AttR^ 



Es^T;i^ + 7^ I dnf + ^^ I dnfiH + i)f). (56) 



V. SHAPE DEPENDENCE OF THE CASIMIR ENERGY, DISCUSSION 



It is instructive to compare the Casimir energy (|56p with the energy in a spherical cavity with equal volume. The 

i?2 



volume and the area of the cavity after deformation are ( up to -^ order ) 



3 ' AttRJ ' 47ri?2 
and 



y=^R'i^~^ fd^f+^J dnf) (57) 



.2,. P f ,r.. . r f .r.s^. 1 



^ = ^"^ (' -tiir^^^tR^ ^"^(' 2^)^) (^^) 



where 

j__a_ . d_ 1 92 



A = — 7^-n0- + -3-^ (59) 



8 

is the Laplace operator on the sphere. The ratio of the energy (|56p and the Casimir energy Eq of the sphere with 
volume ([57)1 is 

^UuoL = i + -r^{fdnfi2 + H)f~^{fdnff). (60) 



This ratio to be examined by its dependence on the shape of the cavity after deformation. One way may be to study 
its dependence on the ratio of the surfaces of the deformed and spherical cavities with equal volumes: 

;||e,. .oi. = 1 - ^ijdnfil + ^A)/ - ^{Jdn/n (61) 

We then express the energy ratio (pO|) as 

^leg. vol. = 1 + ^ J dnf{H - A)/ - 2(|-|e5. ,oL " 1) (62) 

or 

%, jdnfii + iH)f-i,ijdnff s _ 

We know from ([55)1 that the operator H has negative e- values. Thus both of the above relations show that the Casimir 
energy linearly decreases by the increase of the surface. To have a better feeling of this inverse proportionality let us 
consider a simple example. Suppose the deformation function is given by 

fie,(t>)^Piicos9) (64) 

with Pi being the Legendre polynomials. By using ([5^ and ([55]) . we can write ([55)1 for this specific example as 



where the coefficient A is given by 



or by using ([55|) can be written as 



^le,. vol. = 1 - A(0(-^|e9. vol. " 1) (65) 



^<')^^2^WI) '^'" 



^ 2-3,03/3-3,37/2-0,52Z 
^(') = ' 2-1(1 + 1) ■ ^'') 



Thus for large values of I we can write (|65p as 

E~ S 

-rrleq. vol. = 1 — 6/(-— le^. t,oL — 1). (68) 

^0 -JO 

We then can conclude that for less and less smooth deformations we get smaller and smaller Casimir energies. 
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APPENDIX A: THE CALCULATION OF THE VACUUM ENERGY 

Using the uniform expansion formula (|16p for the spectral function (|45p we represent the energy operator (|52p ( 
with jj. = I + 1/2, ui = /ix and e = — ) as 



oo oo J 



H{J) = — V Va^'-'"-'" V Ai,Di{N,J) / dxa;2y2„(a:,£)y2„(a:,£) 



n,m— /—J 

where 



and 



l=J N=-J 

and 



T = 

where 



is the Riemann zeta function 



(Al) 



^"^"'^) = i(tT|^"''"^*^'^^^ ^^'^ 



ti^)- /.,l\l^ , ■ (A3) 



After the change of variables 2n — s + t, 2m — s — t we have 

oo 

H{J)^Y.^'(J) (A4) 

s=0 

where 

^'(■^) = — "E^'"'' y ^-NG-j.ifi) dxx^F,{x,e) 



(A5) 



F,(x,£) = Y. Y,+t{x,e)Y,^t{x,e). (A6) 

By using the Taylor expansion at £ = 

oo „ 

F,(a:,e)=^iyFi")x) (A7) 

k=0 

and the asymptotic expansion (jB13[) we arrive at 

oo 

iJ''(J) =^6^^(2t + 2s-2, J) (A8) 






(AlO) 



br = -± ^Y^^^P r dxX-F^^-^\x) 

and 

(/(A^)) = ^ A5(./(Ar). (All) 



N=-,J 
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For large J we have 

C(2r + 2s-2, J)~ j3-2^-2s (-^^2) 

and 

{Zr+pN-~n - J^^- (A13) 

This asymptotic expressions imply 

H'iJ)-.]^-^', J>1. (A14) 

Thus H'^{J) and H^(J) give the main contributions to the energy operator H{J). Using the summation formulae 
(|BT8)) we get 

HM) = j^{J{J + 1)C(0, J) - ^(^ + y- 1) ^(2, J)) (A15) 



and 



H\J) = ^(ii£±^((2, J) + -^(J + l)4m, J)) (A16) 

^ '^ 2048a^ 2 ^^ ' '^ 16 • 24^ 2^ ^^ ' '^'^ ^ ' 



from which by the virtue of 



1 1 1 

^^^"^^ ~ 7+172 + 2( J + 1/2)2 + 6( J + 1/2)3 + o^-^) 



we read 



9 10 1193 ,,^^, 

^ 1024a 1024a' 1024 • 768a ^ ^ 

APPENDIX B: THE KLEBSCH GORDON COEFFICIENTS 

The spherical harmonics are 

2P^ {cos d) cos mcj), m = l,2,..J 
Y^{9,(b)^{ Piicos9), m = (Bl) 



where Pi{x) is the Legendre function and 



2P;''"'(cose') sin I TO I (/), TO = -1, -2, I 



l-rny.^^ ,, „ d" 



is the associated Legendre function. We use the notation Y^{d,(j)) — Y^{n) where 

n — (cos 6, sin 9 sin 0, sin 9 cos (/)) (B3) 

is the unit vector on the sphere. 

The addition formula for the spherical harmonics is 

Piiiri,n'))= J2 YLin)Yi{n') (B4) 

7n— — l 

where the argument of the Legendre function is the scalar product of the unit vectors n and n'. 
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The Klebsch Gordon coefRcients Kf,,j^ defined by the expansion 



are \ll 



where { fi = 1 + 1/2) 



and 



i+i' 
Pi{x)Pi.{x)^ Y, KiPj{x) (B5) 

J=\i-i'\ 



K^i+n = {J+\WnFi{N,J) (B6) 






Pi (^' J) - .. , N+j -FTr— ^ZT— TWTTTTOT (^^) 



A./ _ ^ ^ V 2 )^\ 2 ) mo\ 



These coefRcients are nonzero if | A^ |< J < 2/ + A^ and iV + J is even number. 
The formula 

n-l 



imphes 



1 B 

lnr(z) = zlnz-z--lnz + lnx/2^+5]^^^^^-^^p^ (B9) 



r(0 + i 1 1 1 , ^ 

In , ,^ ^ -Inz \ BIO 

r(z) 2 82 192z3 ^ ' 



up to the third order in - . The last formula allows as to get the asymptotic expansion of the functions Fi (TV, J) . Up 
to the third order in — for large values of /i we have: 

Fi{N,J)= ^ (i + _L_^) (Bll) 

For the function 

Gi,{^i) = {^i + N)Fl{N,J) (B12) 

we get the following asymptotic expansion 

Gi{^i) = l + — + ^ + ^ + ■■■ (B13) 

where 

2 ' 

J(J + l)-2iV2 

Z3 . Ar ^^'-/J-'"" . (B14) 



Let J is even. Putting J — 2j and N — 2n we have 



^ "^' ^ - (j-«)!(j+n)! 



V^ A J „^Ne _ Y^ i 1 U - » -^ 2 ^^ U -^ " -^ 2^ ^2n9 



Ai^e'^" = }_^ - '■" ,. ^";„; , ^,, ^^ e'""'^ (B15) 

N=-J n=-j 
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In the new variable m ^ n + j we have 






— ' TT m\{2j — to)! 

N=~J m=0 ^ ■' ' 



r2j, 



which is exactly the series representation for zonal spherical functions Y^ ■' (cos 9) 

J 



E A 



5(,e'^» = y/(cos0). (B17) 



N=-J 

The same formula is true for odd J. For example we have 

J 






N=-J 
J 



2 

Af=-J 

■^ 3j3(J + 2) + J(J-2) 



^ ^,^4 ^ o^ w^Y^w-^; ^ ^B18) 



N=-J 
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